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Abstract 

These lectures are designed to introduce the methods and results of large 
N c QCD in a presentation intended for nuclear and particle physicists alike. 
Beginning with definitions and motivations of the approach, we demonstrate 
that all quark and gluon Feynman diagrams are organized into classes based 
on powers of 1/N C . We then show that this result can be translated into def- 
inite statements about mesons and baryons containing arbitrary numbers of 
constituents. In the mesons, numerous well-known phenomenological proper- 
ties follow as immediate consequences of simply counting powers of N c , while 
for the baryons, quantitative large N c analyses of masses and other properties 
are seen to agree with experiment, even when "large" N c is set equal to its 
observed value of 3. Large N c reasoning is also used to explain some simple 
features of nuclear interactions. 



*Lectures presented at "11th Indian-Summer School of Nuclear Physics," sponsored by the Insti- 
tute of Nuclear Physics, Rez and the Faculty of Mathematics and Physics of Charles University, 
Prague, Czech Republic, September 7—11, 1998. 



I. MOTIVATIONS AND FUNDAMENTALS 



The purpose of these lectures is to define, explain, and exhibit the methods of large 
N c QCD for phenomenologists of both the nuclear and particle bent. Several excellent 
reviews on this topic exist in the literature, each with a different emphasis, and to each I 
am indebted for aspects of the pedagogy presented in the current work. Witten's article on 
large N c baryons [Q , one of my favorite papers, actually begins with a review of large N c 
mesons. Chapter 8 in Coleman's Aspects of Symmetry [[| covers the development of large 
N c through 1980. Many of the newer developments with baryons are included in Manohar's 
1997 Les Houches lectures and Jenkins' review [Q. 

I assume that initially many readers may find the notion of large N c to be a rather 
exotic and unlikely approach to studying the strong interactions, as it seemed to me when 
I began to learn it some years ago. As always with an unfamiliar topic, it is important to 
start with clear definitions and motivations. We begin with the assertion that QCD-like 
theories possess a peculiar limit in which the physics of the strong interactions becomes 
much simpler. It is the so-called large N c limit, in which the number of color charges, N c (= 
3 in our universe) is taken as a free parameter (i.e., the gauge group is SU(iV c )), and one 
considers the limit N c — > oo. Physical quantities are then considered in this limit, where 
corrections appear at relative orders 1/N C , etc., the 1/N C expansion. Shortly after its 

discovery by 't Hooft || in 1974, the method was applied to mesons with many qualitative 
successes in explaining phenomenology. However, it has only been in the past several years 
that much success has come in the analysis of baryons, with much of it quantitative. 

However, before a detailed discussion of these very interesting points, it is important to 
consider some simple but incisive questions. How does large N c work? What does it predict? 
And, perhaps most importantly, is it believable? 

In this section, we discuss how and why large N c should work, and how it is implemented 
at the very basic level of QCD, in terms of quarks and gluons. To start with, the very notion 
of the approach creates an interesting paradox: How can increasing the number N c of degrees 
of freedom actually simplify the strong interaction problem? It is certainly true that the 
detailed interactions of any particular quark or gluon become much more complicated if 
one insists on solving the theory exactly. To illustrate this point, let us introduce two 
analogies; we shall come to appreciate that they represent the way large N c works for mesons 
and baryons, respectively. The first analogy is that of a particle moving in a complicated 
potential, one that requires many parameters for a complete description, and the second 
is the many-body problem in mechanics. Both are difficult because of the large number 
of available dynamical degrees of freedom. However, if increasing the number of degrees 
of freedom permits one to consider only collective couplings in a systematic way, then the 
form of the dynamics simplifies. In the first analogy, the potential might be decomposable 
into a simple dominant piece (which nonetheless originates from a large number of more 
fundamental interactions), while the other portions are subleading perturbations. In the 
second, statistical mechanics describes the system's gross features in terms of just a few 
collective quantities, such as temperature or volume; fluctuations about these values are 
highly suppressed. 

With these qualitative arguments in mind, what is special about the large N c limit of 
QCD in particular? 
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1. It is a perturbative QCD expansion. Compare QCD to QED; why does the latter work 
so much better? Obviously, «qed ~ 1/137 is a perturbative expansion parameter, 
while «qcd = 0(1) at hadronic energy scales is clearly not perturbative. At high 
energies, «qcd runs to smaller values due to asymptotic freedom, but one would prefer 
a perturbative parameter that works at all scales. In fact, only one such parameter is 
known to exist for QCD-like theories, namely, 1/N C . 

2. Physics simplifies in the large N c limit. As we show presently, certain classes of Feyn- 
man diagrams dominate over others in physical amplitudes, and this organization leads 
to many interesting simplifications. For example, we shall see that a color singlet qq 
pair always appears at leading order in 1/N C as a single meson. 

3. It provides an explanation of phenomena hard or impossible to explain in field the- 
ory. For example, OZI suppression and the physical dominance of resonant-mediated 
processes are consequences of large N c , as we shall see. 

4. It seems to work, even for N c = 3. This is perhaps the most startling reason; indeed, 
one might expect that 1/N C = 1/3 is too large to be a reasonable expansion parameter. 
However, such an objection may be answered in two ways. One is described by Coleman 
[§] as Witten's "wisecrack": Consider «qed = e 2 /4n ~ 1/137, which leads to e ~ 0.3 ~ 
1/3. If 1/N C is too large an expansion parameter, why don't we say the same thing 
about the QED charge e? Turning the argument around, one can find cases where the 
expansion parameter turns out to be 1/N% ~ 1/10, which is more satisfactory to some 
critics of the approach. However, I prefer a much more pragmatic answer, which is 
just that one should perform a given calculation with N c a free parameter, set it to 3 
at the end of the day, and see for oneself whether or not the factors of 3 are supported 
by the experimental measurements. The ultimate justification of any approach in this 
view is whether it successfully explains physical observations and can be tested in new 
situations. 

Before proceeding, consider the evidence that N c = 3 in our universe; it is important 
to check that there is nothing fundamentally immutable about this particular value that 
would impede consideration of larger N c . Historically, the original reason for the invention 
of color [H] was to explain how fermionic baryons could have a completely symmetric spin- 
flavor-space wavefunction,[] as suggested by experiment. A canonical example is A ++ (J Z = 
+3/2) = (u t)(« T)( M T); which is completely symmetric under spin-flavor indices. Since 
the spatial wavefunction for this low-lying state is expected to be completely symmetric, so 
is the combined spin-flavor-space wavefunction, and therefore one requires a new degree of 
freedom with (at least) 3 values completely antisymmetric under quark exchange. 

The experimental confirmation of the color idea came with the measurement of the ratio 

-Rhadron = cr(e + e~ — > hadrons)/cr(e + e~ — > / u + n~) (1.1) 



Of course, everything said here about flavor is equally true in the two-flavor case of isospin. 
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away from thresholds for new flavor production, which is found, as originally predicted 
0, to be a factor of 3 larger than what one would predict from a colorless quark theory. 
Similarly, perturbative QCD, in those high-energy kinematic regimes where it is believed 
valid, only gives results in agreement with experiment when one sets N c = 3. Furthermore, 
the calculated rate for it — ► 77, which matches experiment, includes a factor of N c = 3. 

Finally, the Standard Model is only renormalizable if the particle content is chosen to 
eliminate the possibility of chiral anomalies. In particular, since quark charges are known 
from deep inelastic scattering experiments to be integral multiples of 1/3, one can show that 
chiral anomalies cancel only for the choice N c = 3. For universes with larger values of N c , 
the cancellation of anomalies is accomplished by adjusting the quark charges accordingly, to 
be integral multiples of 1/N C . 

With these lessons in mind, let us consider the valence structure of hadrons in large 
N c . In order to make progress, we first make the crucial assumption that color confinement 
occurs for arbitrarily large N c ; otherwise, such a universe bears no resemblance to ours, 
since then free colored hadrons or even free quarks can be abundant. Such an assumption 
is not required, of course, since the physical extrapolation from N c = 3 to N c — > 00 might 
lead from a confining to a nonconfining phase at some value of N c > 3, but then the large 
N c limit would be of little phenomenological value since we have no empirical knowledge of 
nonconfining strong interactions. 

Mesons of valence flavor structure 51^2 with color indices a, (3 for arbitrary iV c have the 
normalized wavefunction 



where 5 a p serves to combine the qq pair into a color singlet. The salient point is that mesons 
still have the quantum numbers of a system with a single quark and a single antiquark for 
arbitrary N c . Baryons, on the other hand, are built in analogy with the N c = 3 observation 
that fermionic states with completely symmetric spin-flavor-space wavefunctions require at 
least as many colors as quarks. On the other hand, once we specify that there are N c colors, 
SU(iV c ) group theory requires a minimum of N c colored quarks to form a color singlet. 
Therefore, large N c baryons require exactly N c valence quarks. Given the valence quarks 
qi, q2, ■ ■ ■ , qN c , the normalized wavefunction including color indices a%, 012, • • • , &n c is 



Furthermore, only odd iV c produces fermionic baryons, since each of the N c quarks has spin 
1/2. The large N c limit may be used equally well to consider even-iV c bosonic baryons, but 
of course these are not phenomenologically relevant. 

We are now ready to consider large N c QCD itself, as developed by 't Hooft. To the 
non-initiate who has been told of large N c results but never taught its origins, the whole idea 
might still seem a bit like necromancy. But like any good magic show, the essential predictive 
power of large N c just comes from trap doors (topology) and mirrors (group theory). 

Let us begin with the usual particle content of QCD with the color gauge group SU(iV c ). 
Quarks carry a color index in the fundamental representation of the group, which means 
that the standard representation is a column vector q a with N c components, one for each 







3 



color. Diagrammatically, it is represented by an arrowed line. Antiquarks carry a color index 
in the fundamental conjugate, or anti-fundamental, representation, and therefore appear in 
standard form as a row vector q a with N c components. They are also drawn as arrowed 
lines, but the (anti-) color flow is understood to be opposite the direction of the arrow. 
Gluons, which appear in the adjoint representation of SU(iV c ), carry one color and one 
anti-color index, and appear in standard form as traceless N c x N c matrices A a /3, which 
can be drawn as two parallel lines whose arrows point in opposite directions. This is called 
't Hooft's double-line notation-^ it is actually extraordinarily convenient in capturing the 
color physics, because exact color conservation implies that color lines are continuous, while 
the confinement assumption means that color lines form only closed loops. The three QCD 
vertices in double-line notation are depicted in Fig. p]. 

The double-line notation is particularly useful as a guide to developing an understanding 
of how the QCD gauge parameter, g s (s = strong), scales as N c — > oo. In fact, as 't Hooft 
argued in his original paper ||, QCD Green functions only have smooth, finite limits as 
N c — > oo if g s oc l/y/N~ c . Since this is such a central result in the analysis to come, we 
prove the claim in detail below. However, let us ponder for a moment why such a result is 
interesting. Given that g s oc 1/ ^/W C i one can show that certain classes of Feynman diagrams 
dominate over others, in possessing fewer powers of 1/N C , an organization known as large 
N c power counting. The dominant class of graphs is denoted planar, which means that these 
graphs can be drawn in the double-line notation in the plane of paper, with any quark line 
(if one is present) placed on the perimeter, and with color lines arranged to cross only at 
vertices. Each nonplanarity, as will be shown, costs a suppression factor of l/iV|. Moreover, 
diagrams with no internal quark loops dominate over the others, in that each internal quark 
loop costs a factor of 1/N C . So, in summary, large N c counting provides an organization of 
all Feynman diagrams into a well-defined and countable hierarchy of classes. 

We now present a proof of 't Hooft's scaling of g s . Vertices typically fall into three types. 
There are bilinears that create qq pairs from the vacuum, and are accompanied by g° s \ let 
there be Vi of these current insertions.^ Each of V3 trilinear (qqg or ggg) vertices includes 
gl, while each of V4 quartic vertices includes g 2 s . Confinement requires that the C color lines 
in the digram completely close; then one may choose any of the N c colors to appear on each 
loop, providing a combinatoric factor of . In total, then, the diagram scales as 



Since the diagram has no external lines by color confinement, both ends of each propagator 



2 A couple of notes for the inquisitive. First, double-line notation works as well for Faddeev-Popov 
ghosts (which also live in the adjoint representation), so we may simply think of them in large 
N c counting as another type of gluon. Second, the double-line notation is strictly accurate for the 
gauge group U(iV c ) rather than SU(iV c ). What about the extra U(l) gluon? Since it is only one 
gluon out of N 2 , one can show [|||| that it always gives diagrams subleading by at least two powers 



In addition, the bilinear may also create n gluons, for which gauge invariance requires including 
a factor g™, but such vertices can also be shown to satisfy the counting rules. 



Aoc g 



V 3 +2V A A rC 




s c ' 



of N c . 



4 



P terminate on a vertex: 



2P = Y, nV n = 2V 2 + 3V 3 + W A . 



(1.5) 



n 



Defining the total number of vertices 



v = J £v n = v 2 + v 3 + v 4 



(1.6) 



n 



one has 



Aocg s 



2P-2V N C 



(1.7) 



It is very fruitful to use the double-line notation to interpret the diagram as a polyhedron, 
much like a geodesic roof. Each color loop C becomes a /ace, where each propagator P is an 
edge and each vertex V of the diagram is also a vertex of the solid figure. Then each quark 
loop, which possesses one less color line than if it were replaced by a (double-line) gluon 
loop, may be thought of as a missing face, or topological boundary B, while each nonplanar 
line may be interpreted as a handle (since it represents one edge forced to pass in front of 
another but not meet at a vertex), or topological unit of genus G. We relate these features 
to topological objects so that we may now use a famous theorem due to Euler: For each 
orientableQ figure, there is a topological invariant called the Euler character xe satisfying 



The first term still depends on the particular diagram, while the second has a purely topo- 
logical power. Now observe that, if g s falls off faster than l/^J^S^ c , then diagrams of any 
given topology with the fewest color loops, and hence a minimal number of gluons, dom- 
inate; but this is a trivial theory, in that only valence diagrams are important. Clearly 
the gluon degrees of freedom are vital to understanding the strong interaction, so we reject 
this possibility. On the other hand, if g s falls off slower than 1 / ^/W c , then diagrams of any 
given topology with the most color loops dominate; but this is a nonsense theory, for one 
can always add more color loops to diagrams of a given topology, and in this scenario one 



Double-line diagrams are orientable, i.e., possess a well-defined and continuous definition of 
normal vector to their surfaces, and hence, an inside and outside. This is because each double- 
line edge has one arrow pointing in either direction: Recall that such a feature is used when 
proving Stokes' theorem in calculus. In turn, this feature of the double-line notation comes from 
using SU(iV c ) as a gauge group, where gluons may be thought of as having (quark + antiquark) 
quantum numbers. The same would not be true in, e.g., 0(JV c ). 



Xe = C-P + V = 2-2G-B. 




Therefore, 




i2-2G-B 



(1.9) 
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could make no predictions at all. However, we know that this is not the case, since pertur- 
bative QCD in high-energy regimes works well. We conclude that g s oc l/\fN~ c is the only 
nontrivial, sensible large N c scaling^] giving the result 

A oc N 2 C - 2G ~ B , (1.10) 

meaning that the large N c counting of any particular diagram is purely topological. Since 
this is the case, for any given diagram one may add an arbitrarily large number of planar 
gluons and remain within the same topological class, i.e., the diagram maintains the same 
large N c counting. However, each additional quark loop or nonplanar gluon adds one unit 
to B or G, leading to an additional 1/N C or 1/N 2 suppression, respectively. 

Another (although actually equivalent) way to demonstrate g s oc 1/ \fW c is to start with 
the QCD renormalization group equation, 



. d g 8 (fj) = gsipf 

djj ° 167T 2 



V—, — = -oo T~ a » +0(g s (//)), (1.11) 



where &o = (lliV c — 2Nf)/3 > 0, with Nf the number of light quark flavors, and the 
positivity of b being the reason we know that QCD is asymptotically free. Let us rescale 
g s = g s N~ l l 2+e , with g s finite as N c — > oo. Then ( |1 . 1 1| ) reads 

dg s {n) 1 / N f \ g 3 s (fi) 2e 5 



^ = -i{ ll - 2 t)1^ N '' + 0M (L12) 

We immediately see that, if e < 0, g s is constant and g s instantly runs to zero for N c large 
(assuming that the theory has a perturbative limit), leading to a trivial theory of valence 
diagrams only; on the other hand, if e > 0, g s runs infinitely fast, making a nonsense theory 
with no meaningful predictions. Again, we conclude that only g s oc 1/ ^/W c leads to a 
nontrivial and sensible theory. 

We conclude this section by illustrating in Figs. |2| and |3| the large iV c counting in the 
three different approaches: counting powers of g s and N c , counting features of the equivalent 
polyhedra, and counting topological invariants. These examples should provide the reader 
with hands-on experience with the theorems we have considered. 



II. MESONS IN THE 1/N C EXPANSION 

Having shown that the large N c limit of QCD is only nontrivial and sensible if g s oc 
1 / y/Nc, and that Feynman diagrams involving quarks and gluons fall into classes labeled by 
powers of N c with planar diagrams at leading order, what have we gained? Certainly, the 
observed strongly interacting particles are hadrons rather than quarks and gluons, and so if 
there were no obvious relation between the two pictures — beyond the valence quark model, 



5 Yet again, nature might not require "real" large -/V c QCD to satisfy the conditions of nontriviality 
and sensibility, but if it does not, one cannot extrapolate from N c — > oo back to N c = 3 and hope 
to make predictions. 
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which unlike large N c contains no explicit gluons — large N c would hold little phenomeno- 
logical interest. But in fact one can derive an elegant and simple large N c relationship 
between hadrons and their fundamental QCD constituents. In this section, we consider this 
derivation and its consequences for mesons. 

A good starting point for these discussions is to consider the particle spectrum of large 
N c . In our N c = 3 world, one observes conventional mesons with the quantum numbers of 
qq, baryons with the quantum numbers of qqq, and more recently, hybrid mesons || with 
the quantum numbers of (qq + gluons). On the other hand, one might equally well expect 
to see glueballs, which are defined as mesons containing all gluons and no valence quarks, 
exotic multiquark mesons with the quantum numbers of (qq + extra qq pairs), exotic baryons 
with the quantum numbers of (qqq + extra qq pairs), and hybrid baryons, with the quantum 
numbers of (qqq + glue). Members of the latter class have not yet been observed, and with 
the exception of hybrid baryons, large N c provides an explanation. 

To see how it works, first consider QCD two-point diagrams as depicted in Figs. 
The current J creates a qq pair from the QCD vacuum at some spacetime point, the quarks 
propagate and interact for some interval, and then the pair is annihilated by a current J'. 
We ask, what kinds of hadrons can appear in the intermediate state? Equivalently, using our 
confinement assumption, we want to discover what color-singlet combinations of quarks and 
gluons are revealed if we cut the diagram at any convenient place. We claim that, to leading 
order in 1/N C , cutting such a diagram leads only to a single qq meson — no multiquark 
mesons, no glueballs. 

The proof of this statement relies on planar diagrams with no internal quark loops (which 
alone eliminates multiquark mesons from our discussion) dominating in the 1/N C expansion. 
Typical planar diagrams in the double-line notation are those exhibited in Fig. |2]. The 
unique feature possessed by planar diagrams and no others is that, when cut, each color 
index appears adjacent to its corresponding anticolor index. That is, one cuts through an 
entire color loop before passing to the next one. In terms of fields, the equivalent (perhaps 
nonlocal) operator has color structure 

a A ai A a2 ■ ■ ■ 4 a "-i n an (9 1 ^ 

qa\ A C12 71 a-i A a n q ■ 

In terms of the generators of color indices, one may use the usual Gell-Mann matrices A a to 
expose the underlying structure 

(A ai r a2 (A a2 ) a2 Q3 ---(A a ") Q - 1 Qn . (2.2) 

Now one uses a little SU(iV c ) group theory to complete the proof. A particularly convenient 
way to do this is to exhibit the commutation relations obeyed by the generators: 



\ a ,\ b 



2if abc \ c , {\ a , A 6 } = -^S ab + 2d abc \ c . (2.3) 



Combining these, 

(ATUA b r a3 = ^ ab S a \ 3 + (d abc + if abc ) (A c r Q3 , (2.4) 

where the first term, which is color singlet, is suppressed by 1/N C , while the second term, 
which is color adjoint like the original gluons, is of leading order, 0(N®). By induction, 
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every string of gluon fields is color adjoint to leading order, and therefore the only way to 
make a leading-order color singlet from the diagram is to take the operator qAA ■ ■ ■ Aq as an 
irreducible whole. The conclusion is that only a single meson, and no glueballs or multiquark 
mesons, appear at leading order in 1/N C . 

This important result can now be used to prove a wide variety of consequences. Since 
each planar two-point diagram with a qq pair appears as a single meson, the sum of planar 
two-point diagrams with momentum transfer p 2 must also equal a sum of meson propagators, 

P 

Planar diagram sum = OiNl) = — — - — -, (2.5) 

n P 2 ~ ml 

where f n , the nth meson decay constant, is defined as the transition amplitude between 
meson state n and the vacuum through the current J. Using that p 2 is arbitrary and that 
the equality must hold for any (sufficiently large) value of N c , one finds that m 2 n = 0(N®) and 
fn = 0(y/Nc). Furthermore, for p 2 very large, the diagram must approach its perturbative 
QCD result, ~ ln(p 2 /u 2 ), where u is the renormalization point. However, no finite sum of 
factors fn/ip 2 ~ m n) produces a logarithm in p 2 , so the number of meson states in large N c 
must be infinite. Since one may restrict this calculation to diagrams of a given fixed J PC , it 
is also true that the number of mesons of any particular quantum numbers is also infinite. 

One may continue in this fashion to consider three-, four-, and higher-point functions. 
The diagrams in these cases resemble those of the two-point function, except for the inclusion 
of more current insertions, and thus are all O(N^). However, unlike in the two-point diagram, 
there is now more than one way to collect the quark and gluon fields into color singlets. 
For example, for the three-point function, each of the three current insertions may create 
a meson, which annihilate at a trilinear meson vertex (Fig. |3j). Alternately, one current 
insertion may create two mesons, which propagate to the remaining two insertions. The 
three-point diagram is then given by 

v-^ /1/2/3 • (3-meson vertex) ^ (Amplitude to create 2 mesons) • / 2 / 3 

lib (pi -m i)(P2 -™|)(pI ~ m i) ts (pI ~ m l)(pl ~ m i) 

Using that m 2 n ~ AT?, f n ~ \fN~ci one discovers that the 3-meson vertex scales as l/y/W c , 
while the amplitude to create 2 mesons scales as N®. 

Diagrams with more current insertions allow for more types of internal meson vertices, 
and the N c scaling behavior of each type of vertex can be derived inductively. For example, 
meson vertices from the four-point function can be obtained by using results from the two- 
and three-point functions. In general, given an m-point function, each possible n-point 
diagram, n < m, is required to appear with the large N c scaling derived in a previous 
induction step (a statement of unitarity) with all possible permutations of meson propagators 
(a statement of crossing symmetry). One finds in this way the general results that both the 
m-meson vertex and the amplitude to create m mesons from the vacuum scale as iV c 1_m//2 . 
This important result appears to have been first obtained by Veneziano 0. 

Such a scaling implies immediately that mesons are free, since the propagator, i.e., the 
"two-meson vertex," scales as N®, and therefore neither blows up nor vanishes as iV c — > 00, 
while the 1 / \fW c scaling of the three-meson vertex, and the additional \fW c suppression for 
each additional meson, implies that mesons are stable — a remarkable result! But, could 
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mixing with possible exotic states spoil this reasoning? We now discuss the fate of exotics 
to show that this does not occur. 

First consider glueballs; for simplicity, depict those with only two valence gluons, whose 
two-point diagrams resemble those of Figs. ||]and|| except that the exterior quark is replaced 
by a gluon. Then all the large N c counting is the same as before, except that the single color 
line of the quark is replaced by the double line of the gluon. Therefore, each diagram is a 
factor of N c larger, so that the planar diagrams scale as N%. The same proof as for mesons 
shows that these diagrams are dominated by the one-glueball state, and so one finds that 
the masses and decay constants of glueballs scale as N® and N*, respectively. Continuing 
with the same logic as before, the g-glueball vertex scales as N^~ 9 . 

As for glueball-meson mixing, the simplest planar diagram for the process is exhibited 
in Fig. [|. Converting it to double-line format, we find two color loops and two factors of 
g s , so that the diagram scales as N*. Using our previous results for the glueball and meson 
decay constants of and \/W c , which appear at the current insertions, we see that the 
glueball-meson coupling scales as 1/ y/N~ C i meaning that glueballs and mesons do not mix 
at leading order. Therefore, glueballs and mesons are separately free, stable particles as 
N c — ► oo, with the combined result that the g-glueball, m- meson vertex for m > scales as 

ATl-g-m/2 

Next consider multiquark mesons, e.g., those with valence quantum numbers qqqq. These 
are created by current insertions of the form J = qqqq{x), as depicted in Fig. || To leading 
order in 1/N C , the color loops close either as (1,2), (3,4) or (1,4), (2,3). This means that 
the diagram has two distinct closed color loops and scales as however, we know from 
our previous reasoning that each such loop is a single meson to leading order in 1/N C . To 
create a single meson with qqqq quantum numbers requires an irreducible color loop, such 
as (1,2,3,4) or (1,4,3,2); however, such a diagram scales as and is thus subleading. We 
conclude that multiquark mesons are suppressed in l/N c . 

Finally, we turn our attention to hybrid mesons. From the previous paragraphs, one 
might expect that we shall dispose of these as well. But, in fact, they are present in the 



large N c limit jlOfl . To see this diagrammatically, consider hybrids of valence structure qqg; 



one simple planar diagram representing their propagation is given in Fig. ^a. One finds two 
color loops and two factors of g s , meaning that such diagrams scale as N*, the same as the 
leading conventional valence meson diagrams as in Figs. 0, |3|. In particular, one finds that 
hybrid mesons have masses of 0(N®) and decay constants of 0(y/N^), exactly like those 
of conventional mesons. As for meson-hybrid mixing, one leading diagram is depicted in 
Fig. |7|&, which has two color loops and two factors of g s , thus scaling as AT*. Of course, the 
symmetries of spacetime require that the conventional and hybrid mesons in this case must 
share the same quantum numbers. Using that both conventional and hybrid mesons have 
decay constants scaling as \/~N~ c , one sees that the meson-hybrid coupling is 0(N°), so that 
the two mix freely in large N c . An equivalent observation is that all mesons in large N c are 
expected to have a finite hybrid component, and the large N c limit sheds no light on the 
mystery of why almost all observed mesons have only the quantum numbers of conventional 
mesons. 

It should also be noted that hybrids with exotic (non-gg) quantum numbers (such as the 
newly-observed 1~ + ) do not mix with conventional mesons at leading order in 1/N C , but 
are nonetheless not hard to create. The leading three-point diagrams allow for one or more 
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of the current insertions to have exotic quantum numbers, which means that the process 
of exotic meson production occurs at the same order in N c as that of conventional mesons. 
Glueballs, on the other hand, are harder to create from meson sources and slower to decay 
in the large N c limit; we leave the verification of these facts as an opportunity for the reader 
to experiment with double-line diagrams. 

In summary, our analysis of the large N c meson spectrum teaches that, for the phe- 
nomenology of conventional mesons, one may ignore glueballs, multiquark states, and hy- 
brids with exotic quantum numbers. On the other hand, hybrids with conventional quantum 
numbers should be treated just like conventional valence mesons. The general multiparticle 
vertex for g glueballs, m mesons, and h hybrids scales as 

jyl-g-m/2-h/2 ^ ^.7) 

for m > or h > 0. When m = h = 0, the exponent is 2 — g. The physical mesons, which 
contain a hybrid component, are free, since the propagator scales as N®; stable, since the 

3- meson vertex is 0(l/-y/i\Q, and non-interacting, since meson-meson scattering, i.e. the 

4- meson vertex, is 0(1/N C ). 

In general, one concludes that diagrams with the fewest mesons allowed tend to dominate 
physical processes. This immediately gives a field-theoretic explanation of the dominance of 
resonance-mediated decays and the existence of excited mesons narrow enough to measure. 
Although one could imagine an alternative strong interaction in which every mesonic process 
with a substantial amount of energy release simply produced a spray of pions and nothing 
else, the real strong interaction very definitely produces a rich and complex spectrum of 
observable resonances. Large N c QCD gives an extraordinarily economical explanation for 
this basic observation. 



One success of the large N c limit is the OZI rule ||11|| , that meson processes requiring the 
annihilation or qq pairs, or equivalently, the presence of an all-glue intermediate state, are 
phenomenologically suppressed. The canonical example is the decay of the 0(1020) meson, 
which is believed to have the valence structure ss. It decays into a KK pair, for which there 
is almost no phase space (IrriK ~ 1 GeV) 83% of the time, but only 17% of the time into 
p7r, 7T7T7r, and other nonstrange modes, for which there is much more available phase space. 
Apparently, ss annihilation is suppressed, and large N c offers a very simple explanation: A 
decay like — > KK may be represented with a standard, leading-order (O(N^)) three-point 
diagram as we have discussed, but annihilating the strange quark closes its quark color loop 
and requires another to create the light hadrons in the final state. Since each quark loop 
costs a factor of 1/N C , the leading OZI-suppressed diagrams inherit this suppression.^ 

The observed fact that mesons appear in nonets of flavor SU(3) rather than octets, 
i.e., that the flavor octet and singlet mesons tend to mix and have comparable masses, is 
another success of large N c . The canonical example is 7/(547) and r/'(958), which are the mass 
eigenstates of mixing between a flavor singlet 771 and flavor octet rjs; a priori, the 771 could 
have turned out much heavier than ng. However, one can show that diagrams differentiating 



6 The case with finite N c = 3 is much more subtle, as real OZI suppression appears to have 



interesting dynamical origins above and beyond the mere factor of 1/N C = 1/3: See [12] 
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between their masses either break SU(3), large N c , or both. To see this, begin in the SU(3) 
flavor limit where u, d, and s quarks are not distinguished. Then all planar two-point 
diagrams such as in Fig. contribute equally to both masses, but a diagram with pure glue 
intermediates contributes to m 2 (r]i) alone, since gluons, like rji, are flavor singlets. Because 
the latter diagrams are suppressed by 1/N C , one concludes that terms symmetric in SU(3) 
that distinguish the masses — which a priori could have been large — are in fact suppressed 
by 1/N e effects: 

m 2 ( m ) - m 2 ( Vs ) = O (1/N C , m s ) . (2.8) 

As a final example of mesonic applications of large N c , consider the famous 't Hooft model 
fT3|| . By definition, the 't Hooft model is large N c QCD in 1 space and 1 time dimension 
(1+1). Although it is called a model, it is actually a full-fledged quantum field theory that 
is exactly soluble, in that all hadronic Green functions may be obtained entirely in terms of 
quark degrees of freedom. For example, one can compute meson masses precisely in terms of 
quark masses. Furthermore, the only asymptotically free states in the theory are confined, 
color-singlet hadrons, which appear in an infinite tower of increasing masses. 

So, although we assumed confinement in large N c QCD in 3+1 dimensions, it appears 
to be an immediate consequence in 1+1. Such a result might be expected based upon naive 
dimensional analysis. For, consider the momentum space potential derived from a single 
gluon exchange propagator in the instantaneous approximation: 

V(q) ~ 1/q 2 , (2.9) 

and Fourier transform to position space. In 3+1 dimensions, one has V(r) ~ 1/r, the 
usual Coulomb interaction, which is asymptotically free as r — >• 00; in 1+1, V(r) ~ r, 
which confines as r —>■ 00. That the confining potential survives relativistic and field theory 
corrections is, however, nontrivial. 

The 't Hooft model is soluble precisely because of the two features of its definition. Large 
N c guarantees that nonplanar gluons and virtual quark loops are suppressed, while working 
in 1+1 allows one to choose gauges for the gluon field such that gluon-gluon couplings vanish. 
Such terms arise in minimal substitution from the commutator [A^, A v ] term; working in a 
gauge such that one of the two components of A vanishes guarantees that the commutator 
does also.[] The only remaining diagrams are "rainbow" and "ladder" types (Fig. |8]) and 
these can be summed using "bootstrap" -type Schwinger- Dyson equations, essentially the 
same trick that is used to sum the geometric series. 

Performing this summation for the Green function of a qq system, one obtains the meson 
wavefunction in terms of the 't Hooft equation, 

= (— + T^M tfTW ~ fdy^iy) Pr-^, (2.10) 
\ x 1 — x J Jo (y — x) z 

where the available kinematic invariants are the quark masses M and m, the resultant meson 
mass eigenvalue /1, and the 1+1 analogue x of the deep inelastic variable xbj, defined by 



furthermore, if one chooses a linear gauge, then Faddeev-Popov ghosts also vanish. 
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"mes ' i^mes 

Qualitatively, x and 1 — x represent the fraction of meson momentum carried by the quark 
and antiquark, respectively. The quark masses are renormalized by 

m R = m - —— , (2.12) 

from which we note that the 1+1 strong coupling g s actually has dimensions of mass, and 
using that g s ~ \/\fW C) it is natural to describe masses in units of g s ^jN c /2n. Such a 
quantity acts in the same way in 1+1 as Aqcd in 3+1, in that it distinguishes "heavy" from 
"light" quarks. The 't Hooft equation is known to possess an infinite tower of solutions, 
labeled by n — 0, 1, 2, . . ., which alternate in parity (P = (— l) n+1 ). As n —>■ oo, the meson 
masses scale as fi^ — > n 2 n. Since x is bounded between and 1, the 't Hooft equation is 
qualitatively similar to the quantum mechanical problem of a particle in a box. 

The utility of the 't Hooft model is that it can be used to study difficult problems of 
QCD that have not been fully solved in 3+1, such as the meson form factors, deep inelastic 
scattering, and the saturation of hadronic rates by partonic diagrams. For example, consider 



T4| the nonleptonic weak decay of a meson containing a heavy quark into lighter mesons. 
As the heavy quark mass M is increased, the light antiquark and gluon degrees of freedom 
should exert less and less influence on the decay of the heavy quark, and eventually the 
process should be well described by the diagram of free heavy quark weak decay. Since one 
can compute both the hadronic and partonic decay rates in the 't Hooft model, the two 
can be compared directly. The result is Fig. |9], in which one sees that the two rates rapidly 
approach one another for large M. 



III. BARYONS IN THE 1/N C EXPANSION 

Recalling from Sec. I that baryons for quarks of N c colors have the quantum numbers 
of N c quarks, one is immediately faced with the dilemma of how to understand the physics 
of a system with N c —>■ oo constituents. A particularly simple and elegant paradigm for 
dealing with this situation is the mean- field Hartree-Fock baryon picture due to Witten [|T| , 
in which each quark moves in the field caused by the other N c — 1 quarks acting collectively. 
To lowest order in 1/N C , this field is static — after all, a system of N c — 1 quarks possesses 
a great deal of inertia compared to a single quark. Then each quark in the ground state 
satisfies the same field-theoretical wave equation (which, for very heavy quarks, reduces to 
the Schrodinger equation); although the exact form of this wave equation depends upon the 
details of QCD and is therefore unknown, much can be said from its existence and its large 
N c scaling properties. 

For starters, since each quark in the ground state satisfies the same wave equation, each 
one has the same wavefunction 4>(r), which we claim scales as N®. For the scaling claim to be 
sensible, the potential due to the other N c — 1 quarks in which each of the N c quarks moves 
must also scale as N®, or equivalently, baryon diagrams must never give interaction energies 
more than O(N^). To see that this is true, one must first understand large N c counting for 
baryon diagrams. The potential comes from diagrams with connected interactions such as 
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those depicted in Fig. [LC]. In addition to factors of N c from g s oc 1/ \fN~ c and color loops, 
one has combinatoric factors due to the presence of N c quark lines. For example, diagrams 
involving two quarks can be chosen in v c C*2 = N C (N C — l)/2 = O(N^) ways, diagrams 
involving three quarks can be chosen in n c C% = N c (N c — l)(N c — 2)/6 = 0(N%) ways, and so 
forth. The three diagrams in Fig. [10] are all seen to appear at O(N^), and so are among the 
leading-order diagrams. It is not hard to prove that no diagrams are of lower order, and that 
internal quark loops and nonplanar gluons (See Appendix |A]) cost factors of 1/N C , meaning 
that a consistent large N c counting scheme compatible with the Hartree-Fock picture does 
indeed exist for the baryons. 

Since <p(r) is the same for each quark in the baryon, it follows that each quark has 
the same charge distribution and occupies the same space; therefore, baryons remain the 
same size as N c — *■ oo. In contrast, large A nuclei grow in size with A because of the 
phenomenological saturation of nuclear forces. Whereas large N c baryons and large A nuclei 
both have masses proportional to the numbers of constituents N c and A, respectively, large 
N c baryons are like rigid containers with radius oc N® and density oc N^, while large A 
nuclei are like close-packed spheres with density oc A and radius oc A 1 / 3 . This is why the 
Hartree-Fock picture works for baryons but not nuclei. 

Within the Hartree-Fock picture it is possible to obtain large N c behavior for common 
physical processes. For example, consider (Fig. [12]) baryon-baryon (BB) and baryon-meson 
(BM) scattering. In the BB case, one may choose from among N c quarks in each baryon, 
while at least one gluon is typically required to keep both baryons color singlets and conserve 
momentum, providing a factor (g s / ^/W c ) 2 . The amplitude for the process thus scales as N*. 
In the BM case, the counting proceeds as before, except that one loses one of the iV* 
combinatoric factors, meaning that the amplitude scales as N®. Do these factors make 
sense? For BB, since baryon masses are O(N^), an amplitude of O(N^) leads, via the usual 
textbook formulae, to a nonvanishing, finite (0(N®)) cross-section. For BM, since meson 
masses are 0(N®), the cross-section formulae lead to the scattering of the meson — but 
not the baryon — with 0(N®) probability; the physical picture in this case is the meson 
scattering from a fixed potential center.^] 

Although the Hartree-Fock picture provides a simple archetype for the large N c baryon, 
its value is limited since it only describes the leading-order N c counting for a given process. 
This was also the case for the mesons, but now with N c constituents in the baryon, it 
becomes possible to perform studies of subleading effects in 1/N C , which involve some subset 
of the quarks acting collectively. In the next two subsections, we shall discuss the two known 
techniques for studying the subleading 1/N C structure of baryons: consistency condition and 
spin-flavor algebra methods. Both make use of an operator basis acting upon the baryon 
states, and so the bulk of the 1/N C suppressions is carried by the operators. The remainder of 
this lecture is dedicated to explaining these two schemes and their phenomenological results. 
The consistency condition approach turns out to be related closely to an operator analysis 



To be a bit more precise about the calculation, these results are most clear if one uses the covariant 
normalization for baryons of E/M particles per unit volume. Although a mere bookkeeping device, 
this factor is 0(N®) in the large N c limit (i.e., not infinite like 2E) and therefore gives the standard 
result for scattering from a potential center in the BM case as N c — > oo. 
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based on the Skyrme model |15| , while the spin-flavor method turns out to be closely related 
to operator analysis performed in the quark model. 



A. Consistency Conditions 

The starting point of this scheme is to consider baryon-meson scattering not in terms of 
quarks, but rather hadrons. Suppose that the meson is a tt, which, as a Nambu-Goldstone 
boson of chiral symmetry, couples to baryons derivatively through the isovector axial vector 
current A^ a = 57^75 r a q, where r is an isospin Pauli matrix. The scattering process consists 
of two insertions of A^ a on the baryon line, which can occur in two orderings as depicted in 
Fig. |13[ Since each quark in the baryon can couple to A^, one expects its matrix elements 
to be 0(^)11] and the full vertex is written as 



B' 



qrWg 



B) d^/U = ~gN c X™, B d t Ti a /U, (3.1) 



where the coupling constant g and the operator X m , each of 0(N®), are defined by this 
equation. Note that only the space components of A^ a appear in the leading term of the 
nonrelativistic expansion of the quark bilinear, and only these are needed here, since at 
leading order in 1/N C the baryon may be taken at rest throughout the process. This implies 
also that the ir scatters elastically from the baryon, and so in the first diagram in Fig. [13] the 
baryon is off-shell by an amount E n , whereas in the second the amount is (— E n ). Therefore, 
since the derivatives produce factors of pion momenta k, the amplitude for the two diagrams 
together is 



M 



iN 2 c ~g 2 k l k'3 

J 7T 



B'B 



(3.2) 



Using that f n oc \fW C) M. naively appears to scale as N^, rather than N® as concluded 
previously. One evades this conundrum by noting that many baryon states may appear 
on the intermediate line, i.e., are inserted between the two X's in the commutator; then 
cancellations can render the matrix elements of the commutator 0(1/N C ) even though matrix 
elements of each X individually are 0(N®), as pointed out by Gervais and Sakita in 1984 
||16|| . Such cancellations are called consistency conditions. 

The mathematics behind this approach was developed in the 1960s as one of many 
methods to study strong interactions, and is generically called the method of induced repre- 
sentations of a contracted spin-flavor algebra. In order to explain this intimidating term, we 
first consider what is meant by a contracted algebra. In essence, it means that one or more 
generators of an algebra is scaled so as to change the commutation relations. For example, 
perform & l/N c expansion of the operator X: 

X ia = X ia + ±_ x ia + . . . ? (3 3) 



9 The "seagull" diagram, where the two pions meet at one point on the baryon, is suppressed since 
it couples to the isospin quantum number of the baryon, which is taken to be 0(N®) for nucleons. 
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in terms of which we have seen that 



Via -yjb 
A > A 



0. 



(3.4) 



By virtue of its defining indices % and a, X m is a spin-1, isospin-1 tensor, and therefore by 
definition satisfies 



J l ,X 



■ ijk ykb 
it A , 



Also, as usual, 



J\J J 



ie ijkjk 



ja rb 



i\xi b 



ie abcjc 



ie abc Xl c . 



r,j l 



o. 



(3.5) 



(3.6) 



which are the defining relations of spin-flavor SU(2) x SU(2). How do the commutation 
relations with X alter the algebra? To see what happens, compare the algebra defined by 
(pP-(p|) with that of spin-flavor SU(4). In addition to containing SU(2) x SU(2), SU(4) 
also has commutation relations with the combined spin-flavor generator G ia = J 1 <S> I a , 



J\G ja 
G ia G jb 



ie ijk G ka , 



/ a , G j 



ie abc G ic 



- ( e ijk 5 ab J k + e abc 5 ij I c ) . 



Note that the two algebras become the same if one writes 

X l Q a = lim G ia /N c , 



(3.7) 



(3.8) 



and works to leading order in 1/N C , so that the right hand side of the commutator [G, G] 
vanishes. This is the contraction of the spin-flavor algebra. One immediately concludes 
that, due to the contraction, the smaller spin-flavor symmetry SU(2) x SU(2) is promoted 
to the larger symmetry SU(4), at leading order in 1/N C . For three light flavors u, d, s, for 
which the flavor symmetry is SU(3), one finds SU(2) x SU(3) — > SU(6) at leading order. 
This enlargement of the symmetry due to the contraction allows one to find a number of 
new relations among physical quantities. In fact, the SU(6) spin-flavor symmetry jI7| has 
been known to work well for baryons — but rather poorly for mesons — since the 1960s; in 
large N c , one sees at last a field-theoretic reason how such a symmetry might arise. 

Next, a representation of an algebra is induced when one attempts to specify the quantum 
numbers of a given state as classified according to a generator of the algebra, and discovers 
that other quantum numbers must be specified to label the state completely. This happens, 
for example, when one classifies the irreducible representations of the Lorentz group; once a 
given inertial frame (i.e., set of boost coordinates) is chosen, one finds that there still exists 
a "little group" of transformations needed to specify the representation, which correspond 
physically to rotations in that frame: Choosing the frame induces the full representation. In 
the present case, once one specifies eigenvalues of the operator X , one finds 



18 that there 



remains a "little group" of generators with quantum numbers K } k, like those of angular 
momentum, K = 0, 1, 2, . . ., and k = —K, —K + 1, . . . , K — 1, K. That is, choosing Xq 
induces an irreducible representation specified by the states | X , K,k). 

Of course, X is not an intuitively familiar physical operator, so one then projects the 
states | X , K, k ) onto the space of states | J, J3, J, J^,K ). Certainly isospin and spin are 
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familiar, but what of the induced quantum number Kl As shown by Dashen, Jenkins, and 



Manohar ||19|| , K turns out to have the same properties as the number of strange quarks in 
the baryon, K = N s /2. The relation mentioned earlier between the consistency condition 
approach and the Skyrme model, namely, that the group theory associated with the Skyrme 
model is precisely that obtained from the consistency condition method, is also discussed in 
that work. 

The consistency condition ( |3.4| ) and others like it derived from multiple pion or kaon 
scattering with nucleons provide fertile ground for phenomenological analysis. One simply 
plugs an operator possibly contributing to a given physical process into the consistency 
relations, sees if it is satisfied, and if not, discards the operator. This process builds a 1/N C 
expansion for a given operator, and hence for physical quantities. For example, one can 
derive the pion-baryon coupling constant relation 



9nNN/9nNA = 3/2 + 0(1/N%), 



(3.9) 



or show that X[ a in (|3.3|) is proportional to Xq q , meaning that effectively 



X 1 



0(1/N ( 



2 ) 

C 11 



(3.10) 



so that nonstrange axial currents have no 0(1/N C ) corrections, or show that the masses of 
nonstrange baryons of spin J in the same spin-flavor multiplet, such as the nucleons and 
A's, are first split at relative order J 2 /N^ 



01, i.e., 



M = N r 



(3.11) 



where cq and c\ are unknown coefficients of 0(N®). These relations arise simply be- 
cause the operators that would violate them are disallowed by the consistency condi- 
tions. Phenomenology seems to agree well with this analysis; in the examples given, 
9-kNa/ 9-kNN = 1-48 experimentally, while the relative mass difference between the nu- 
cleon and A is about (1232 - 940)MeV/|(1232 + 940)MeV = 0.27, which is very close 
to [J A (J A + 1) - Jn(Jn + l)]/iV c 2 = 0.33 for N c = 3. 

A large body of additional work has been performed using the consistency condition 



approach. To sample but a few, baryon magnetic moments []2T|, orbitally excited p2| , |23" 
and hybrid p3 have been studied in this scheme. 



B. Spin-flavor Algebra 

In the consistency condition approach just discussed, we saw that a full combined spin 
and flavor symmetry arises from the separate spin and flavor symmetries to leading order 
in 1/N C , due to the group contraction. Why not, then, describe all operators acting upon 
the baryons in terms of operators with definite SU(2F) properties, where F is the number 
of light quark flavors, acting upon the sum of the quarks? That is, the operator basis in 
SU(6), for example, is defined as 
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T a = qi(t ® y) q a , 



G m = d (j ® yJ q a , (3.12) 

where a 1 are the usual Pauli spin matrices, A a denotes Gell-Mann flavor matrices, and the 
index a sums over all N c quark lines in the baryons. One can readily show that J, T, and 
G satisfy an SU(6) algebra analogous to the SU(4) algebra Eq. (|3.7| ) -H 

An important point to realize is that, although baryons are described in this approach 



in terms of N c "quarks," a large N c constituent quark model |25| is not assumed. Indeed, 
recalling the lesson from Sec. |lj that valence and hybrid mesons mix freely, one might well 
expect the same for hybrid baryons |[24|| . However, the only assumption really used to 
perform this operator analysis is that baryons for large N c continue to have the same total 
quantum numbers as those obtained from an N c quark system. Then, a "quark" in the 



sense of Eq. (|3.12 ) is not the same as the dynamical entity of the QCD Lagrangian, but 



rather a group-theoretical object representing one part in N c of the physical baryon. The 
two definitions coincide when the quarks are taken very heavy, so that their interactions 
with the internal gluons become less important in determining the baryon structure. 



It was shown by Manohar |26| that results of the Skyrme model obtained from its operator 
structure (and hence obtained from the consistency condition approach) and those obtained 
from the operator structure of the quark model (and hence obtained from the spin-flavor 
algebra approach) coincide as N c — > oo. Therefore, either approach is equally good from the 
mathematical point of view for describing large N c baryons. Moreover, the total result for 
any given quantity must be the same regardless of which approach one uses, which is simply 
the statement that the physical result must be independent of mathematics. One sees that 
the consistency condition and spin-flavor algebra methods differ for a given physical quantity 
only by a reorganization of 1/N C corrections. 

The initial work on the baryon 1/N C expansion in the spin-flavor algebra approach was 
performed by Carone, Georgi, and Osofsky |27 ] and Luty and March- Russell [^8| . The basic 



results of these analyses is that many of the old SU(6) static relations between baryonic 
matrix elements are recovered by means of an operator expansion in 1/N C . That is, any 
given physical operator that scales as N™ may be written in the form 

o^^TTj^r^r^y-', (3,4, 



^ ~"\N J \N ) V TV 

n,p,q \ iv c/ \iv c / \ iV c 



To be precise, one replaces I a — > T a , e — * f , and the [G, G] commutator becomes 

[gta = l _ 5 ijjabc T c + l_ 5 ab e ijkjk + l_ e ijk d abc G kc_ (3^3) 
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where it is understood that spin and flavor indices are to be contracted in such a way as to 
agree with the transformation properties of 0^ m \ 

The factor of 1/N C accompanying each J, T, and G reflects the fact that the matrix 
elements of these operators can often add coherently over the N c quarks and thus can be as 
large as O(N^). In fact, this was the central problem with the spin-flavor algebra method: 
One hopes to describe a given quantity with some small number of operators, but if a 
number of the operators in (|3.14|) have matrix elements of O(N^), where can it justifiably 



be truncated? Before addressing this problem directly, it is useful to consider the nature of 
the baryon multiplets in this approach. 

Since we are making use of representations of the symmetry group SU(2F) for opera- 
tors acting upon the baryons, the most convenient description for the baryons is in terms 
of multiplets of SU(2F). Recall that this is the symmetry of independent rotations in 
all spin SU(2) and flavor SU(F) coordinates, under which the ground state baryons ap- 
pear to fill a completely symmetric representation. It is very convenient to describe the 
representations in terms of Young tableaux, in which each quark, which appears in a fun- 
damental representation of SU(2F), is represented by a single box. Boxes arranged in a 
row indicate symmetrization of the corresponding SU(2F) indices, while boxes arranged in 
a column indicate antisymmetrization. The tableau for the completely symmetric SU(2F) 
multiplet of N c quarks is a row of N c boxes, as in the first line of Fig. IT3]; using standard 



group theory techniques, this representation can be decomposed into a tower of component 
SU(F)fl avor x SU(2) spin representations, as indicated in the succeeding lines of Fig. [14]. Each 
Young tableau in this tower has N c boxes. If you are familiar with the old SU(6) theory of 
baryons, then you will recognize the special case of N c = 3, F = 3 as the decomposition 

56= (8, J = -) © flO, J= -V (3.15) 



2J V 2 

The significance of this result is twofold: First, the number of SU(F) x SU(2) representations 
increases as N c is increased from its physical value of 3, and second, the individual SU(F) 
multiplets themselves become much larger as N c is increased beyond three.Q This point is 



illustrated in Figs. ^ and ^ for the case of three light flavors; note how the familiar SU(3) 
octet (which contains the nucleons) and the decuplet (which contains the A resonances) are 
obtained from shrinking the representation down to iV c = 3. 

Let us now return to the conundrum of where to truncate the expansion ( |3.14j ). We 
see that working with N c large generates an enormous number of baryon states, but we are 
only interested in the phenomenology of those with which we are familiar, i.e., that already 
exist for N c = 3, which we denote physical baryons. Where do they appear within the 



large multiplets of Figs. 15, 16? The mathematics alone allows for a great deal of freedom 



for placing the physical states, but the usefulness of ( E>.14| ) as a perturbative expansion in 



1/N C imposes constraints on this choice. For example, since the expansion breaks down if 
J = 0(Nl) ox I = 0{Nl), we take physical baryons to have J = 0(N°) and / = 0(iV c °). As 



There is one exception, namely, F = 2. In that case, all of the two-box columns become 
isosinglets and may be ignored, and then the isomultiplet representations are independent of N c 
— indeed, they satisfy I = J. 
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for the strangeness content, consider for example the top-right entry of the multiplet Fig. [TJ|; 
its spin-flavor structure is uud in a total spin-1/2, isospin-1/2 state (just as for an N c = 3 
proton) plus (iV c — 3)/2 ud pairs, each combined into spin-0, isospin-0 quantum numbers. 
It is tempting to identify this state as the N c > 3 analogue of the proton, and one may 
do this by fiat; indeed, in light of our statements above, the physical proton should have 
J = I = 1/2 even for N c > 3. Apart from the top-right state, the physical proton might 
occupy any of a number of sites in Fig. [T5] that lie directly below this one. However, each 
of these states contains a number of s quarks, which are heavier than u and d quarks, and 
therefore these baryons should be heavier than those on the top row. If the physical proton 
contains some valence s quarks, it decays weakly to some other lighter baryon, a process 
that opens up a phenomenological Pandora's box. One does not encounter any of these or 



similar problems if one identifies the physical baryons to lie at the tops of Figs. [15| and |1C . 

Even taking these assignments into account, some of the operators, such as T 4 ' 5,6 ' 7 and 
G m , can still give matrix elements of 0(N C ), even when evaluated only on the physical 
baryons. All is not lost, however, since we have not to this point made full use of the 
mathematics. Since the operators can be thought of as a set of vectors spanning the space 
of physical observables, it must be that there are only as many independent operators as 
independent observables. One can obtain a number of operator reduction rules based 
upon both the structure of the spin-flavor algebra SU(2F) and properties particular to the 
completely symmetric ground-state representation. 

To see how this works, collect the SU(2F) generators into a uniformly normalized set: 

A A = { (V /2 <g> l) /VF, (11 ® X a /2) /y/2, V2 (0-72 ® A a /2) } -> Tr A A A B = -5 AB . (3.16) 

Then two simple operator reduction relations are 

ql (11 ® TL)q a = N c t, (3.17) 

(<zi aV) (?} aV) = + 2F ) i 1 ~ a - ( 3 - 18 ) 

Equation ( [3. 171) is actually nothing more than the total quark number operator in the baryon, 
which of course produces a factor N c . The left-hand side of Eq. ( |3.18| ) by definition is the 



quadratic Casimir of SU(2F) in the same way that J 2 is the Casimir of spin SU(2), and 
therefore is also proportional to the identity operator. Note in both cases how each reduction 
by one operator leads to an additional power of N c \ it is this feature that maintains the 
consistency of large N c counting within the operator reduction rules. Once the complete 
set of such rules is known, one may produce a minimal set of independent operators to any 
given order in the 1/N C expansion to describe a given physical quantity. 



Analysis of the mass spectrum [30] provides a good illustration of these ideas. Let us 
consider here for simplicity only the eight isospin-averaged masses N, E, 5, A, £*, H*, Q 
of the ground-state spin-flavor multiplet, where the baryon label stands for its mass. For 



N c > 3 the relevant baryon multiplets contain many additional states, as a glance at Figs. 15 



16] verifies, but for sake of phenomenology only these 8 physical baryons need be considered. 
One must therefore generate the leading 8 independent operators in the 1/N C expansion, 
taking into account all operator reduction rules, in order to obtain a set of mass operators 
that forms a complete basis (in the sense of a vector space) for spanning all possible values 
of masses. The operators with attendant 1/N C factors are 
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n c i, r 8 , i-j 2 , -L re-, i- (t 8 ) 2 , j 2 t 8 , ^rcv, JL (t 8 ) 3 . (3.19) 

The mass Hamiltonian operator is simply a linear combination of these operators, each with 
an arbitrary coefficient. One may now compute the matrix element of each operator for each 
baryon, obtaining in this way an expression for each coefficient in terms of baryon masses. 
For example, the coefficient of (T 8 ) 2 is fixed by the combination 

35(2iV - E - 3A + 22) - 4(4A - 5£* - 22* + 3ft). (3.20) 

Our mass Hamiltonian has not yet taken into account all the known physics of the ground- 
state multiplet; in particular, the masses of the baryons are known to break along the 
hypercharge (strangeness) axis, which corresponds to the explicit flavor = 8 indices appearing 
in ( |3.19| ). This explicit SU(3) breaking is known to be relatively small, at the order of e ~ 
0.25-0.30, since the relative mass difference between, say, iV(940) and £(1190) is roughly 
this size. Every time a flavor = 8 index appears, let us include a factor of e in the mass 
Hamiltonian. 

One would now like to compare the matrix element e 2 ((T 8 ) 2 ) /N c to the mass combi- 
nation ( |3.20| ), but both of these expressions still possess scale independence. That is, both 



( |3.20| ) and the coefficients of the mass Hamiltonian may be multiplied by arbitrary fixed 
numbers. One overcomes this ambiguity by comparing to the O(N^) common mass of the 
baryons. On the operator side in our example, this means that the relevant quantity is 
e 2 /iV 2 , which for e = 0.25 and iV c = 3 is about 0.69%. On the mass side, one organizes 
a given linear combination of masses into the form LHS = RHS, where all numerical co- 
efficients on either side are positive, and since the combination vanishes in the spin-flavor 
symmetry limit, the total of the LHS and RHS coefficients are equal. One then forms the 
combination (LHS + RHS)/2. In our example, this is 

- [35(2iV + £ + 3A + 22) + 4(4A + 5£* + 22* + 3ft)] . (3.21) 

The desired scale-independent combination is then |LHS — RHS|/|(LHS+RHS). Here, this 
is experimentally 0.37 ± 0.01%; thus the unknown coefficient of the operator e 2 (T 8 ) 2 /iV c 
has the very natural size 0.53. Had we ignored the factors of N c , the coefficient would 
have been 0.06, its small size suggesting that some important physical suppression had been 



ignored. Carrying out this analysis for all 8 masses gives the plot in Fig. [L7]; here one 
sees graphically that not only the factors e but also the factors of N c = 3 are needed to 
understand satisfactorily the mass spectrum. 

The truly interesting feature of this analysis is that large N c is used to obtain the relevant 
operator expansion represented by ( |3.19| ), but the phenomenological analysis then employs 
N c = 3. One concludes that, at least for this case, the methods of large N c seem to survive 
extrapolation all the way to the rather small observed value of N c = 3. Similar analysis 
in the quark operator representation has been performed for the magnetic moments |31|j32]| 
and axial- vector couplings P2"| , |3"3"|] , where evidence for the predictions of the 1/N C expansion 



is found. 

However, this is not to say that the 1/N C expansion is perfect! An interesting example 
of a case where the expansion for N c = 3 does not explain everything is that of nonstrange 
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1 excited baryons fl34| |, such as the AT(1535). An analysis similar to that described 



above (somewhat more complicated due to the additional presence of the orbital angular 
momentum operator t) shows that a number of the mass Hamiltonian coefficients have 
central fit values that are quite small, albeit with substantial uncertainties. It is not yet clear 
whether this is due entirely to large experimental uncertainties in the mass determination 
of these resonances,^ or additional dynamics beyond the 1/N C expansion, or both. On 
the other hand, if any of the coefficients had turned out too large, it would have been an 
unmitigated failure of large N c , since in that case the claimed 1/N C suppressions simply 
would not be supported by experiment. 



C. Nuclear Physics 

Large N c nuclear theory is a subfield still in its infancy, but the few results obtained 
thus far are encouraging. Interestingly, all of the current results tend to follow from the 
observation that G ia has O(N^) matrix elements for the physical baryons, while those of I a 
and J 1 are 0(N®). In Ref. P5fl , the dominance of G m is used to show that the 1/N C leading- 
order central nucleon-nucleon (NN) interaction is SU(4) symmetric. To see this, consider 
the collected SU(4) generators A A , as defined in Eq. ( p.lb] ). The leading SU(4) symmetric 
operator for the system consisting of nucleons 1 and 2 is 



A- 



(3.22) 



which contains 



G l 



G' 



(3.23) 



where the 0(1/N^) corrections come from and J\-Ji terms. Since G ta is a leading-order 
operator, one concludes that the leading-order interaction is SU(4) symmetric. This reason- 
ing is very similar to that of Sec. |III A| , where group contraction promotes separate SU(2) spin 
and SU(2) isospin symmetries to SU(4), since it depends upon terms in [G, G] (Eq. ( |3.13| )) 
being subleading in l/N c . 

One very interesting consequence of this result is an explanation for the observed Wigner 
supermultiplet symmetry. Many nuclear interactions seem to obey a symmetry among the 
states in the multiplet 



(p T, Pi, n t, n |) 



(3.24) 



and large N c provides an explanation |^| : Given the NN Lagrangian only up to dimension-6 
operators, one has 



C 6 = -\c s (iVtiv) 2 - ~C T (iVV'JV 



(3.25) 



However, data from CLAS at Jefferson Lab should improve upon the experimental values. 
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where the Wigner symmetry is preserved by the first term and broken by the second. How- 
ever, the second breaks SU(4) and therefore is subleading in 1/N C . 

The leading 1/N C behavior of the non-central terms is derived in Ref. 
nonrelativistic, elastic NN potential is expressed through the decomposition 



36 , where the 



NN 



where 



V ° + V>! • <x 2 + V® S L • S + V°S l2 + V°Q 12 
V* + Vl a- x • <x 2 + Vis L • S + V* S 1 2 + V% Q i 2 



Tl • T 2 , 



(3.26) 



S12 = 3 crx • f a 2 ■ f — crx • cr 2 . 
Qu = ^ j>i • L, <j 2 • L} . 



(3.27) 



Again using the dominance of G ta , one can show that the potential terms V °, V}, and V? 
appear at 0(Nl), while every other term is relatively suppressed by at least two powers of 
1/N C . The immediate consequence for the NN interaction is that these terms should domi- 
nate the phase shifts, and indeed this is the case. As has been long known, the conspicuous 
features of the NN interaction are a spin-singlet central potential, for which V^ is prominent, 
and an isotriplet Lorentz tensor force, for which and are prominent. In phenomeno- 
logically successful one-meson exchange potential models these features are represented by 
including large coupling constants for n, a, and uj mesons for the central force, and 7r and p 
for the tensor force. However, the salient point is that the relative importance of these terms 
in the potential derives directly from large iV c counting, regardless of the physical origin of 
these couplings. Any phenomenologically successful model must obey the same pattern. 



D. Further Advances 

Phenomenological studies of the 1/N C expansion comprise an increasingly large volume 
of literature. While these notes are designed as lectures rather than a complete review, I 
would like at least to mention a number of additional interesting problems that have been 
tackled to date. 

• That the 1/N C expansion and chiral perturbation theory of baryons may be combined 
was demonstrated by Jenkins [0. The feature discussed earlier that mesons appear 
in nonets rather than octets of SU(3) flavor (i.e., that the r( meson should be included 
in the chiral Lagrangian) is supported by this construction. One finds in comput- 
ing loops with baryon lines in this Lagrangian (which requires including towers of 
intermediate states: See, e.g., Fig. |13D that the tree-level counting of powers of 1/N C 
is maintained; that is, loops produce complicated functions of the chiral symmetry 
breaking parameters only (such as m^lnm^), not N c . 

• The famous Heavy Quark Effective Theory |38| (HQET) relates physical amplitudes 
of hadrons with heavy quarks (b or c) of different spins or flavors based on the idea 
that a sufficiently heavy quark is like a static color source, the cloud of light quarks 
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and gluons surrounding it barely noticing its particular spin or flavor quantum num- 
bers. The relation of different amplitudes, in particular form factors, means that an 
approximate symmetry is at work in HQET; however, since large N c induces an ap- 
proximate symmetry as well (spin-flavor SU(2F)), one would expect that large N c and 
HQET together produce additional constraints. This is indeed the shown by 

Chow [j39| : There are relations between the form factors for the semileptonic decays 
A fe -> A c £u, S b -> S c £z/, and -»■ £*£z/. 

Large N c and HQET can be combined into a larger and highly predictive effective 
theory, as shown by Jenkins [ilOfl . The large N c spin- flavor symmetry in this case 
is SU(2Fheav y ) x SU(2Fi ig ht) -Q One obtains from this study a hierarchy of charmed 



baryon masses, in analogy to that of |30| , whose masses agree with experimental values. 
One can then use the charmed masses to predict the yet unmeasured beauty baryon 
masses to high accuracy. 



The decays of £ = 1 baryons to ground-state baryons through pions |^2J and through 
photons jR| have been investigated in the 1/N C expansion. These decays will be 



measured in great numbers at Jefferson Lab. 

This sampling of different directions of research should demonstrate that the field of 
large N c phenomenology is quite vital and holds the promise of producing many further 
interesting, relevant results. I invite you to add to the sum of this knowledge. 



IV. CONCLUSIONS 

Large N c QCD is an elegant theoretical construct that organizes all QCD Feynman 
diagrams into a countable set of classes based on their topological properties. The 1/N C 
expansion thus obtained is the only known way to turn QCD into a truly perturbative theory 
at all energies. Predictions of large N c are obtained by counting the numbers of powers of N c 
associated with each diagram; since this hierarchy persists when infinite numbers of Feynman 
diagrams are summed to give diagrams for hadronic processes, the approach graduates from 
one of mathematical curiosity to real phenomenological relevance. 

For mesons, many phenomenologically observed but poorly understood properties simply 
follow from the leading-order term of the 1/N C expansion, while the operator expansion for 
baryons provides a vehicle for systematic quantitative analysis of masses, couplings, and 
other physical observables. This is an active area of current research at the time of this 
writing. 

All of the results presented here, with the exception of the 't Hooft model in Sec. [TT], rely 
solely on the large N c counting rules. Essentially, the counting rules provide nothing more 
than an organizing principle for contributions to physical quantities, while the complicated 
details of the dynamics reside in the coefficients of these operators. 
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For comparison, the .Fheavy = 1 case was explored in [41]. 
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Considering the numerous successes of large N c , it seems altogether possible that the 
eventual solution of QCD will depend intimately on this remarkable limit. This is an opinion 
shared by many of those who perform much more formal field-theoretical work and are 
endeavoring to compute the dynamical coefficient in from of the leading term in 1/N C - 
that is, to sum all the planar graphs. Can this be done for real QCD? If so, perhaps buried 
in the result of this important calculation are treasures like the origins of confinement and 
chiral symmetry breaking. 
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APPENDIX A: "PLANARITY" IN BARYON DIAGRAMS 

The question of what constitutes a "nonplanar" gluon for baryons is more tricky for 
baryons than mesons. Consider Fig. [11]; in Fig. |TT| a one sees an ordinary "ladder" diagram 



that for the mesons would be planar, whereas the crossed diagram Fig. |Tlp would be nonpla- 
nar. For baryons, however, the double-line notation shows that a has no closed color loops, 
whereas b has one. Therefore, since both diagrams share the combinatoric factor n c Ci and 
four factors of g s , the two diagrams scale respectively as N® and iV*; Figure |Tl"| & is therefore 
a leading-order diagram. Indeed, we leave it as an exercise for the reader to draw the color 
flow for the two diagrams, in order to be convinced that maintaining the direction of the 
arrows on the color lines requires the double lines of gluons to be twisted. Since the two 
diagrams differ by the crossing of a gluon line, one may call a "nonplanar" and b "planar" . 

Why has this complication occurred? The key is that all of the quarks in the baryon have 
the same direction of color flow. For mesons, quarks and antiquarks (which have opposite 
directions of color flow) always appear in pairs, and therefore adjacent color loops in the 
double-line notation always have lines pointing in opposite directions. This is precisely what 
is needed to produce orientable surfaces (see Footnote [|) and a convenient diagrammatic 
topological classification of meson diagrams, and what is lacking in the baryon case. Whereas 
for mesons each nonplanar ity cost a factor of 1/N%, we see from our example that nonplanar 
l/iV* suppressions for baryons are possible. 
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FIGURES 



nV 



nV 



FIG. 1. The three types of QCD vertex in standard Feynman notation in the top row, accom- 
panied by their representations in double-line notation in the bottom row. 




FIG. 2. Examples of large iV c counting for Feynman diagrams. Each row gives three equivalent 
representations of the same diagram, the first in standard form, the second in double-line notation, 
and the third as a topological figure. In terms of the quantities in Eqs. fll.7| ) and (|1.§|), the first 
diagram has g° s N^, C = 1, P = 2, V = 2, G = 0, B = 1; the second has g 2 s N^, C = 2, P = 5, 
V = 4, G = 0, B = 1; the third has ^iV c 3 , C = 3, P = 8, V = 6, G = 0, B = 1. All three scale as 

and are planar. Crosses indicate current insertions connecting quarks to the QCD vacuum. 
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FIG. 3. More examples of large N c counting for Feynman diagrams. The first line has g^N^, 
C = 2, P = 8, V = 6, G = 0, B = 2 and scales as iV c °; the second has g 6 s N*, C = 4, P = 10, 
V = 7, G = 0, B = 1 and scales as TV*; the third has c^ 1 , C = 1, P = 8, F = 6, G = 1, B = 1, 
and scales as iV" 1 . Note that only the second diagram is planar and hence leading in 1/N C . 
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FIG. 4. The two distinct types of meson structures equivalent to the QCD three-point function. 




FIG. 5. One leading diagram for glueball- meson mixing. 



28 




(a) (b) 

FIG. 7. Leading two-point diagrams relevant to (a) free hybrid meson propagation, and 

hybrid-conventional meson mixing. 




FIG. 8. Nontrivial "rainbow" and "ladder" diagrams summed by the 't Hooft model. 
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L , , 1 , , , , 1 , , , 

5 10 
M 

FIG. 9. Full rate for nonleptonic weak decay of meson containing a heavy quark of mass M 
as a function of M. The solid line is computed in the 't Hooft model and sums over all allowed 
exclusive modes, while the dashed line is due only to the free heavy quark weak decay diagram. 
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FIG. 10. Three sample leading-order connected baryon diagrams contributing to interaction 
energy. Each square bracket indicates the presence of N c quark lines. Counting combinatoric 
factors, color loops, and powers of g s oc \/N~ c , one finds that each diagram scales as TV*. 
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FIG. 11. Simple ladder (a) and crossed-ladder (b) gluon exchange diagrams for baryons. 
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FIG. 12. Sample 1/N C leading-order diagrams for baryon-baryon and baryon-meson scattering, 
respectively. Brackets indicate quarks collected into distinct hadrons. 



\7c(k) rc(k')/ X^OO n(F),s 



B B' B B' 

FIG. 13. Baryon-pion scattering diagrams used to obtain the amplitude (|3.2[) , and ultimately, 
the consistency condition Q3.4|). 



J = \ 



j — He 

J ~ 2 



FIG. 14. The completely symmetric SU(2F) iV c -box Young tableau, corresponding to ground- 
state baryons, decomposed into SU(F)fl avor x SU(2) sp i n representations. 
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1 1 

12 1 
12 2 1 
1 2 2 2 1 
1 2 2 2 2 1 
1 2 2 2 2 2 1 
12222221 
122222221 
1222222221 
111111111 

FIG. 15. Weight diagram corresponding to the ground-state spin-1/2 representation of SU(3) 
for N c > 3. Indicated are the numbers of states at a given site. The longer side has (N c + l)/2 
sites. 

1111 

1 2 2 2 1 
1 2 3 3 2 1 
1 2 3 4 3 2 1 
12344321 
123444321 
1234444321 
1233 3 3 321 
12222221 
1111111 

FIG. 16. Weight diagram corresponding to the ground-state spin-3/2 representation of SU(3) 
for N c > 3. Indicated are the numbers of states at a given site. The longer side has (N c — l)/2 
sites. 
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FIG. 17. Experimental values of isospin-averaged mass combinations relative to the common 
baryon mass, as described in the text (a ratio called experimental accuracy here). Each one corre- 
sponds to a particular operator in the 1/N C expansion, whose 1/N C and SU(3)-breaking suppressions 
e label each point. 
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